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The effect of disorder on magnonic transport in low-dimensional magnetic materials is studied 
in the framework of a classical spin model. Numerical investigations give insight into scattering 
properties of the systems and show the existence of Anderson localization in ID and weak localization 
in 2D, potentially affecting the functionality of magnonic devices. 


The propagation of spin waves [T] is in the focus 
of modern research because of its importance for spin 
caloritronic applications mg and for future information 
processing devices which might either rely on magnonic 
[7] instead of electronic transport or combine electronic 
with magnonic transport. The dynamics of spin waves 
is mostly described by the Landau-Lifshitz-Gilbert equa¬ 
tion, a nonlinear equation of motion that describes the 
wave propagation as well as some degree of dissipation, 
included phenomenologically either following Landau and 
Lifshitz [H] or Gilbert [g. This dissipation limits the co¬ 
herent wave propagation to a spatial scale set by the 
propagation length C, that depends on the material prop¬ 
erties, especially the damping constant [ini E]. The 
microscopic origin of the damping is inelastic scattering 
with, e. g., phonons m- Static imperfections of the mag¬ 
netic crystal, on the other hand, induce elastic scattering, 
which has two effects. First, it turns ballistic into diffu¬ 
sive transport, and, second, it might suppress transport 
completely, as first shown by Anderson in 1958 for spin 
diffusion in disordered lattices m- 

Meanwhile it has been established for many different 
kinds of waves that quenched disorder in combination 
with phase coherence can lead to a complete suppression 
of transport, confining eigenmodes to spatial regions of 
a finite extent given by the localization length ^loc [la¬ 
in addition, there is also the weak-localization regime 
where diffusive transport still prevails, but mesoscopic 
effects of phase-coherent scattering can be observed m- 
Arguably, the most famous phenomenon of that kind is 
coherent backscattering (GBS), an effect that relies on 
long-range phase coherence and can therefore be seen as 
a gauge of the microscopic processes that eventually en¬ 
tail Anderson localization. GBS produces an enhanced 
intensity for the elastic scattering of an excitation with 
wave vector ko into the opposite direction —kg, and has 
been directly observed with, e.g., light HMg , acous¬ 
tic [El E], seismic [E]? as well as matter waves |E] . 
In contrast, localization phenomena for spin waves have 
been studied rather scarcely, mostly in amorphous mate¬ 
rials with random anisotropy piHE] and by analogy with 
hard-core boson excitations on disordered lattices |E 1 - 
Eg. But since different types of defects are very common 
in magnetic crystals, it is important to study their conse¬ 


quences for magnonic transport. Indeed, a localization- 
induced breakdown of regular transport would severely 
hamper the functionality of real-world devices. It is the 
purpose of this paper to study localization phenomena 
with spin waves on the basis of numerical calculations. In 
particular it is important to determine the length scale 
of Anderson localization, and compare it with dissipative 
mechanisms which also limit magnonic transport unmi]. 
Localization effects are known to be most relevant in low 
dimensional systems. We therefore investigate strong 
localization in one dimension and GBS in two dimen¬ 
sions. Also, experimental setups include often thin films 
or nanowires, which might be treated as low dimensional 
materials as soon as their extent is of the order of or lower 
than the wave length of the magnons. 

We study transport of magnons modeled as classical 
spin waves within the framework of an atomistic spin 
model m- The model comprises dimensionless mag¬ 
netic moments (called “spins”) = li^/ps on a lattice 
of sites r*, I = 1,... ,N, with the local magnetic mo¬ 
ment and /is the reference value for the magnetic mo¬ 
ment in the clean material. Each spin S* has neigh¬ 
bors S'", m = l,...,A^nb at positions a™ relative to site 
1. The Hamiltonian of the system realizes a Heisenberg- 
type spin model 
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with ferromagnetic exchange interaction J > 0, a uniax¬ 
ial anisotropy constant and an external magnetic field 
B*. £; = 0,1 states the occupation of a site, which will 
be needed for defects. The spins evolve in time according 
to a torque equation. 
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where 7 is the magnitude of the gyromagnetic ratio. This 
equation of motion corresponds to the Landau-Lifshitz 
(LL) equation in the limit of vanishing damping, where 
the total energy is conserved. The microscopic time scale 
of this model is tj = ps/^J (« 50fs for iron). The 
natural order of magnitude for distances is the lattice 
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FIG. 1. Strong localization of a spin wave packet, launched with initial wave vector fco = 0.3/a (in units of lattice spacing a) 
and initial width cro = 50/\/2 o in a disordered ID chain with a local field B present on a fraction q of random sites, (a) RMS 
wave packet spread (Jr = [(r^) — as function of time for a defect field strength of B = 0.4 J/^s- At long times, magnonic 

transport comes to a halt, (b) Asymptotic in-plane magnetization profiles show exponential localization over the localization 
length ^loc. 


constant a ~ lA, and for local magnetic fields J/^ts « 
100 T. 

In the following we take the system to be globally mag¬ 
netized along z by choosing a small anisotropy 0 < <C 
J. Small-amplitude linear excitations, known as spin 
waves or magnons, are then confined to the xy-plane, 
s’’ = Sl — iSy, such that the local wave intensity 
measures the in-plane magnetization. It will also prove 
fruitful to analyze the momentum-space density |5kp5 
where 5k denotes the Fourier transform of SK In the 
clean, simple cubic, d-dimensional lattice under consid¬ 
eration, the magnon dispersion reads 


Wk 
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For infinitesimal anisotropy dz J, these spin waves are 
the gapless Goldstone modes of the ferromagnetic phase. 

The discrete translation symmetry of the lattice is bro¬ 
ken by the presence of defects. We consider two kinds of 
uncorrelated defects, distributed with number density g 
on randomly chosen sites rA One kind is a local mag¬ 
netic field B-^ = (0,0, B) along the easy axis, with B* = 0 
everywhere else. (So here Si = 1 VZ.) The other kind is 
non-magnetic substitutional disorder with missing mag¬ 
netic moments, Sj = 0 at defects sites and with ei = 1 
everywhere else. Since we investigate finite systems, cal¬ 
culated quantities usually depend on the defect config¬ 
uration, and all results presented below will include an 
ensemble average, noted (...). 

We study the out-of-equilibrium, long-time spin wave 
dynamics, by integrating the equations of motion Q nu¬ 
merically, using an implicit Adams scheme. The initial 
condition at time t = 0 is a quasi-monochromatic Gaus¬ 
sian wave packet 

5o = Aexp [iko-r'- (r'- ro)^/4CTo] (4) 


with amplitude A (in the range 0.01... 0.1) and rms 
width (To around the initial position Tq, launched with 
finite wave vector kg into the bulk disordered lattice. 

In a first step, we study a one dimensional spin chain, 
where disorder should manifest as strong localization. 
We model defects by a random local field along z, taking 
a finite value B with probability g and vanishing with 
the complementary probability 1 — g. (While the other 
defect model of local missing spins is arguably more re¬ 
alistic in the bulk, we do not consider it for the ID case 
since it breaks the exchange coupling and thus trivially 
confines the excitations to disconnected segments.) Our 
ID simulation results are summarized in Fig. Initially, 
the rms width of the wave-packet spreads in time, as 
shown in panel a). At longer times, the width saturates, 
and the spreading comes to a complete halt, which is a 
hallmark of localization. As a rule, the higher the de¬ 
fect density, the stronger the localization effect, and the 
smaller the final extent. The localization scenario is fur¬ 
ther corroborated by the asymptotic in-plane magneti¬ 
zation profiles, plotted in panel b) for different values of 
defect concentration and defect strength on a log-linear 
scale. All profiles show the same characteristic expo¬ 
nential decrease, and a fit to the expected asymptotic 
form exp(—|r — rol/d^ioc) yields the localization 

length ^loc. The conventional factor 4 here emphasizes 
that the ensemble-averaged intensity decays more slowly 
than the typical (i.e. most probable) intensity, which is 
approximately log-normal distributed and decreases as 
exp(—|r — rol/^ioc) [SH Note that the localization 
lengths found in our simulations are at the order of 10^- 
10 ^ a (« 0.01-0.1pm) for the chosen parameters and are 
therefore far below magnon propagation lengths lim¬ 
ited by dissipation, which can be in the range of 10^ a 
(« 1pm) for small damping constants a = 10“^ ... 10“^ 

m- 

The height of the wings depends on the distribution 
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FIG. 2. Clean dispersion relation a;k, Eq. ([^, and elastic 
scattering rate I/cjicTs for a particular set of wave vectors 
inside the first Brillouin zone of a simple cubic lattice. Data 
points are obtained for a missing-spin defect density of = 
0.1 and are averaged over 20 disorder configurations. The 
scattering rate vanishes at the central F-point. Strong 
scattering occurs for intermediate wave vectors. 


details near the center, which are found to deviate from 
the simple exponential cusp predicted in Ref. [33] . This 
is not surprising, given that our simulations do not match 
the assumptions of the analytical calculations. Notably, 
the wave packet starts with finite initial velocity and cov¬ 
ers a certain range of momenta and energies. Moreover, 
the disorder parameters situate the simulation far from 
the perturbative regime. In particular, ^loc cannot be 
expected to be given by the lowest-order term of an ex¬ 
pansion in the defect strength of independent scatter- 
ers. Interestingly, deviations from the profile predicted 
by Ref. [33] have also been observed in numerical simu¬ 
lations of matter waves in uncorrelated on-site disorder 
[55] . In any case, the strong disorder prevents the sys¬ 
tem from reaching its equilibrium configuration (all spins 
aligned along z). Instead, an in-plane magnetization re¬ 
mains forever written into the spin chain, thus highlight¬ 
ing the lack of ergodicity, one of the chief manifestations 
of localization |57H5D] 

In a second step, we analyze magnon scattering in a 2D 
disordered lattice with the aim of assessing weak local¬ 
ization effects. Since in magnetic materials non-magnetic 
defects are rather common, we place zero spins on ran¬ 
domly chosen lattice sites. The only parameter describ¬ 
ing the degree of disorder is therefore the defect density, 
or percentage of defect sites, which we take to be p = 0.1 
for the data presented below. In order to gain a better 
understanding of the microscopic scattering processes at 
work for this type of disorder, we study the magnon dy¬ 
namics in k-space by evaluating /k(t) = (|5k(t)P). The 
initial wave packet, Eq. Q, is a peak of width aQ^ cen¬ 
tered at ko. Due to elastic scattering off the defects, the 
initial wave packet is depleted, and the peak height de¬ 
creases as exp(—t/Ts), where Tg is the elastic scattering 
mean free time. This characteristic time can be measured 
by a fit to the observed exponential decay, thus revealing 
whether scattering can be qualified weak (wk^s ^ 1) or 


must be considered strong (wkTs ~ 1). 

Fig. plots the dispersion, Eq. ([^, together with the 
reduced scattering rate l/wkTs for selected wave vectors 
in the first Brillouin zone. Near the symmetry point T in 
the band center, the scattering amplitude from k to k' is 
proportional to k • k', as characteristic of p-wave scatter¬ 
ing. According to Fermi’s Golden Rule m, considering 
that the density of states in 2D is constant at low en¬ 
ergy, one should expect the scattering rate to vanish like 
|k|^. And indeed, a quadratic behavior of l/wkTg around 
the origin is consistent with the data. For comparison, 
we also determined the scattering rate for the local mag¬ 
netic field defects, which produce an isotropic, s-wave 
scattering amplitude. Consequently, the scattering rate 
decreases more slowly with |k|, but vanishes nonetheless. 
At first sight, this is at odds with Ref. [21] when trans¬ 
lated from 3D to our 2D case with its constant density 
of states. But, as k —>• 0, the independent-scatterer ap¬ 
proximation used in Ref. m breaks down, and collective 
scattering from impurity clusters eventually leads to a 
vanishing scattering rate, as expected from general prin¬ 
ciples for Goldstone modes at low energy [Hj. The local 
maxima of the scattering rate roughly halfway through 
the band, signaling strong scattering, may be traced back 
to a van Hove divergence of the density of states at fre¬ 
quency u) = shifted down in energy and rounded 

by the disorder. 

For further investigation of the strong-scattering re¬ 
gime, we take cto = 150/-\/2a and ako = (0,0.567r) where 
Tg « 1.59 tj (see highlighted point in Fig. |^, such that 
the elastic scattering mean-free path 4 = |vk| Tg « 3.1 a 
(vk = 9kWk is the group velocity). Due to multiple elas¬ 
tic scattering, partial wave amplitudes appear in modes 
k with the same frequency Wk = Wko, up to a disorder¬ 
broadening of order Thus, we can follow the pro¬ 

gressive, diffusive redistribution of wave vectors over the 
energy shell. For a classical random-walk model, i.e., 
phase-incoherent propagation, one would expect a homo¬ 
geneous distribution over all accessible modes, as a con¬ 
sequence of ergodicity. A failure of ergodicity, instead, 
should show up as distinctive features in the wave-vector 
distribution Ik- 

Fig. i a) shows Ik{t) at a time t = 10 tj ^ Tg well 
in the diffusive regime. Above the diffusive background 
of height /bg, it clearly features the GBS peak at —ko, 
whose presence proves that the memory of the initial 
condition is preserved for very long times. Starting 
from a pure plane-wave excitation, the peak contrast 
C = (/-ko ~ /bg)//bg with respect to the background 
should be exactly unity and constant in time. This sig¬ 
nal has to be convolved with the initial k-space distribu¬ 
tion following from Eq. (|^, and therefore is expected to 
decrease as [12] 


m = 


4 ( 7 ^ 


3a2 


’it)" 


(5) 




4 





t/tj 


FIG. 3. Spin-wave intensity 7k = in 2D k-space (average over 800 defect configurations) at times (a) t = lOtj and 

(b) t = 1000 tj. The initial wave packet can still be seen in (a) as the narrow peak at ako = (0, O.SBtt). The CBS peak at —ko is 
well visible above the diffusive background, distributed along the energy shell Wk = • The width and contrast C of the CBS 

peak decrease in time, (c) Time evolution of the contrast for different settings, evaluating the impact of a weak nonlinearity 
{A = 0.1) and finite damping (a > 0). The dashed line shows the diffusive prediction, Eq. ([^, for the linear case. Whereas the 
damping does not affect the CBS contrast, nonlinearities induce dephasing and reduce the CBS contrast noticeably. 


where the diffusive spread cr^(t) = crg-|-|vk|^ increases 
linearly in time. In the case investigated here, we deter¬ 
mine the transport time Ttr using a Green-Kubo relation 

(Vk • Vko) (t )Vk • Vko oc exp(-t/Ttr), (6) 

'^kGB.Z. 

where Af = (time-independent) normaliza¬ 

tion. For the chosen parameters we obtain Ttr ~ 1.37j. 
Fig. c) shows the observed contrast, together with the 
prediction (§, for different simulation parameters. The 
agreement is excellent for a small spin wave amplitude 
A = 0.01. For a larger amplitude A = 0.1 the CBS con¬ 
trast decays faster, a behavior that we attribute to the 
dephasing caused by the nonlinearity (such effects are 
known for, e.g., light m and matter waves m)- 

The third set of parameters includes finite damping. 
Indeed, more realistic dynamic models for spin waves in¬ 
clude the Gilbert damping via an additional term in the 
Landau-Lifshitz-Gilbert (LLG) equation, proportional to 
the damping constant a [3T]. The observed contrast of 
the CBS peak, for the chosen value of a = 1 x 10“"^, re¬ 
mains unchanged compared to the undamped case. This 
is in agreement with the reciprocity principle, well known 
in optics m, namely that uniform damping lowers the 
overall intensity, but preserves the CBS contrast com¬ 
pared to the background. 

In summary, we have numerically studied the influ¬ 
ence of random defects on the propagation of classical 
spin waves in low-dimensional disordered magnets. We 
find evidence for strong (Anderson) localization of spin 
waves in a one dimensional spin chain. In a two dimen¬ 
sional disordered lattice, a clear coherent backscattering 
signal proves the presence of weak localization effects— 
a well-known precursor for Anderson localization. These 
findings underpin the importance of defect-related effects 
on magnonic transport, and define new limits for the 


propagation of spin waves in addition to the usually as¬ 
sumed Gilbert damping. In dimensions higher than one, 
the crossover to the strongly localized regime is hard to 
reach by direct numerical integration because it typically 
occurs at much longer times and for much larger system 
sizes. If one tries to increase the fraction g of missing 
spins too much, the lattice becomes disconnected at the 
percolation threshold, and the excitations become triv¬ 
ially confined to the percolation clusters, which was not 
the regime of interest here. In order to compute local¬ 
ization lengths in the linearized regime (together with 
critical properties of possible localization-delocalization 
transitions in higher dimensions), the method of choice 
is a transfer-matrix approach combined with a finite-size 
scaling analysis m HZ]. Nonlinearities, on the other 
hand, generically suppress the onset of Anderson local¬ 
ization, and lead to subdiffusive behavior instead |48j . 
The quantitative investigation of such effects in substitu- 
tionally disordered magnets poses interesting challenges 
for future work. 
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